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Abstract 


The  expression  u » c r/  ' ^ (where  c is  a constant) 
sometimes  assumed  for  the  displacement  around  a point  imperfec- 
tion (interstitial  or  substitutional  impurity,  lattice  vacancy) 
gives  a non-zero  stress  at  the  surface  of  the  solid.  The  addi- 
tional ’image  displacement1  necessary  to  insure  that  this  stress 
vanishes  is  usually  neglected,  but  may  be  important.  For  example, 
it  accounts  for  from  30  to  50  per  cent  of  the  volume  change  pro- 
duced by  such  defects.  This  and  other  effects  of  the  image  term 
are  discussed.  Miller  and  Russel  have  pointed  out  that  a point 
imperfection  near  the  center  of  a sphere  the  apparent  volume 
change  deduced  from  measurements  of  the  x-ray  lattice  constant 
is  greater  than  the  geometrical  volume  change.  It  is  shown  that 
the  reverse  is  true  when  the  defect  is  near  the  surface,  and 
that  for  a large  number  of  defects  scattered  uniformly  through 
the  sphere  the  geometrical  and  x-ray  expansions  are  equal.  This 
result  is  true  with  or  without  the  image  term,  but  to  establish 
the  equality  of  geometrical  and  x-ray  expansions  for  a body  of 
arbitrary  shape  the  image  terms  must  be  included. 

The  contributions  of  the  image  term  to  various  phenomena 
are  different  and  so  they  cannot  be  merely  absorbed  in  the  'strength 
c of  the  imperfection. 


Introduction 


As  a model  for  the  distortion  of  a crystal  lattice 
by  a substitutional  or  interstitial  atom  or  a vacant  lattice 
site  a number  of  authors  have  taken  a center  of  dilatation  in 
an  isotropic  elastic  continuum.  In  many  cases  calculations 
have  been  made  using  the  elastic  displacement  appropriates  to 
a center  of  dilatation  in  an  infinite  medium.  This  displace- 
ment cannot  be  correct  since  it  would  give  a non -vanishing 
stress  at  the  free  surface  of  the  body.  It  is  convenient  to 
regard  the  actual  displacement  as  the  sum  of  the  displaef  ment 
in  an  infinite  medium  and  an  ‘image*  displacement  due  to  the 
perturbing  effect  of  the  boundary. 

It  is  the  object  of  this  paper  to  show  that  the  image 
term  cannot  be  neglected  in  certain  applications.  1'n  particular, 
if  it  is  omitted,  the  volume  change  due  to  one  or  more  point  singu- 
larities will  be  underestimated  by  a factor  of  1.5  for  metals  or 
1.8  for  alkali  halides.  Again,  if  the  effects  of  the  free  boun- 
dary are  ignored,  any  non-spherical  body  uniformly  filled  with 
point  singularities  would  suffer  a change  of  shape  as  well  as  a 
change  of  size.  When  the  boundary  effect  is  taken  into  account 
there  la  a change  of  volume  without  change  of  shape.  Moreover, 
for  this  case,  a calculation  of  the  change  of  x-ray  lattice  con- 
stant gives  a result  consistent  with  the  geometrical  deformation 
only  if  the  image  displacements  are  included.  Thus,  unusually,  the 
boundary  effects  simplify  the  problem  instead  of  complicrting  It. 

2 . A Point  Singularity  in  an  Infinite  Body 

The  elastic  displacement  due  to  a center  of  tii.<  :< at  ,„n 


in  an  infinite  medium  i.si 


^ A.  H,  Love,  Mathematical  Theory  of  Elasticity.  (Cambridge  University 
Press . 1924). 


vrhere  c is  a constant,  the  ’strength*  of  the  singularity.  {'Je  use 
the  affix  so  to  emphasize  that  (1}  is  only  valid  in  an  infinite 
medium).  The  displacement  has  the  same  form  as  the  field  of  an 
electrostatic  point  charge.  The  real  justification  for  taking  (1)  as 
a rougjh  description  cf  the  elastic  field  of  a point  imperfection  in 
an  infinite  crystal  which  has  been  idealised  as  a homogeneous  isotropic 
elastic  continuum  is  that  it  is  the  only  spherically  symmetrical 
displacement  which  satisfies  the  equations  of  elasticity  and  does 
not  increase  with  . Still,  it  is  convenient  to  have  some  sort 
of  detailed  elastic  model.  As  a model  of  an  interstitial  or  sub- 
stitutional atom  we  might  take  an  elastic  sphere  of  radius  (1  •*  £ 
forced,  into  a spherical  hole  of  radius  :ln  an  infinite  block  of 

the  same  material.  It  is  easy  to  show  that  for  'f~c.  the 

displacement  is  given  by  (1)  with  c " dT'0(  1 -*  0 ' ) /3(1  - (f  ) and 
that  for  nT  < T'  there  is  a uniform  compress  Ion,  { c/  is  Peis- 

son*s  ratio).  The  surface  of  the  hole  suffers  an  outward  displace- 
3-- 

ment  c/Y~0  » increasing  the  volume  within  it  by  4*b* 

This  model  must  not  be  taken  too  literally.  The  lat- 
tice con stunt  of  gold  is  increased  by  a little  dissolved  silver.  The 
misfitting  sphere  mc-dei  would  suggest  that  therefore  a little 


3 


dissolved  gold  would  decrease  the  lattice  constant  of  silver, 
and  this  is  not  true*  All  that  we  can  hope  for  is  that  a par- 
ticular type  of  singularity  in  a particular  matrix  will  be 
characterized  by  a constant  c which  can  be  used  consistently  to 
describe  various  phenomena. 

Again,  a foreign  atom  in  a lattice  exerts  forces  on  its 
neighbors  differing  from  the  * standard  * forces  they  would  exper- 
ience in  the  perfect  lattice.  When  we  assimilate  the  lattice  to 
a continuum  the  standard  forces  are  absorbed  in  the  elastic  pro- 
perties of  the  medium,  but  the  additional  forces  due  to  the  im- 
perfection are  still  outstanding.  Hence  as  an  alternative  to  the 
sphere-in -hole  model  we  may  take  a cluster  of  point  forces  to  re- 
present a point  lattice  imperfection  in  the  elastic  approximation 
For  a sufficiency  symmetrical  relation  of  the  interstitial  or  sub 
stitutional  atom  and  the  lattice,  the  forces  will  be  equivalent 
to  three  equal  ’double  forces  without  moment*  at  right  angles. 

The  displacement  will  then  be  of  the  form  (1)  which  can  be  consi- 
dered1 as  due  to  a distribution  of  body  force 

IM  £ ) - - 12*cKgjrad  v ( f ^ (2) 

(K  is  the  bulk  modulus). 

In  a less  symmetrical  case  (e.g.  carbon  in  iron)  the  cluster  of 
forces  will  be  equivalent  to  three  unequal  double  forces  without 
moment  and  the  resulting  displacement  will  be  more  complicated 
than  (1).  However,  it  is  probably  not  justifiable  to  introduce 
this  refinement  without  at  the  same  time  considering  the  aniso- 
tropy of  the  material,  and  we  shall  not  consider  it  further. 


t 

4 

The  displacement  (1)  can  be  produced  by  heating  a point 
of  the  elastic  medium,  assumed  to  be  non-conducting.  This  can  be 
seen  by  imagining  that  the  misfit  in  the  first  model  was  produced 
by  heating  a sphere  with  originally  no  misfit,  or  in  terms  of  the 
second  model  by  noticing  that  the  thermal  stress2  due  to  a temper- 
ature distribution  T { / ) is  the  same  as  the  stress  produced  by 

2 

S»  Timoshenko,  Theory  of  Elasticity,  McGraw-Hill  Book  Go, , Inc,, 
Nevx  York,  N,  I,  I9T4T, 


a density  of  body  force  proportional  to  grad  T,  which  for  the  hot 

<s 

spot  T » const  '"  ) would  be  of  the  form  (2),  This  analogy  will 
be  useful  in  discussing  a body  containing  a large  number  of  imperfec- 
tions. 

In  general  the  stress  associated  with  a displacement  u^ 
whose  Cartesian  components  are  u^,  U2,  is 

pij  " ^ \ j div  “ * n (3) 

\ i / 

where  X and  p.  are  Lame's  constants.  From  the  electrostatic  inter- 
pretation it  is  at  once  clear  that  the  divergence  and  curl  of  (1) 
vanish,  so  that  the  stress  produced  by  the  singularity  when  in  an 
infinite  body  is  simply 


3» 


pij 


(4) 


A Point  Singularity  in  a Finite  Body  with  a Stress-Free  Surface. 


Consider  now  a center  of  dilatation  in  a finite  body  with 
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a free  surface*  First  mark  out  the  surface  S of  the  proposed  body 
in  an  infinite  block  of  material  and  introduce  the  singularity. 

The  displacement  is  given  correctly  by  (1).  Across  any  surface 
element  of  S there  is  a stress*  p^n^dS,  where  n^  is  the  normal  to 
S.  Thus  if  we  remove  the  material  outside  S the  displacement  will 

* 

Throughout  the  paper  we  use  the  convention  that  a repeat es  suffix 
is  to  be  summed  over  the  values  1,  2,  3* 

continue  to  be  u*  only  if  we  apply  a distribution  of  surface  trac- 

at 

tion  pijnj  to  S,  Removing  this  distribution  to  give  a body  with  a 
stress-free  surface  is  equivalent  to  applying  an  additional  distri- 
bution The  elastic  state  will  then  be  given  by 

H.  * 'C  * X.1*  "ii  ' * Plj  <5> 

where  the  image  stress  p^  is  the  stress  which  surface  tractions 

-pjjnj  would  produce  in  the  body  and  the  image  displacement  is 
is  related  to  pjj  by  (3).  pjj  is  clearly  free  of  singularities 
within  S and  satisfies 

(Pij  ♦ ) nj  « 0 on  S.  (6) 

T 

Unlike  u , ux  has  in  general  a non-vanishing  divergence. 

SK  'VW 

If  the  surface  traction  and  surface  displacement  are  pre- 
scribed {they  are  not,  of  course,  independent)  the  elastic  field 
Inside  the  body  can  be  found  by  integration.  When  only  the  surface 
traction  is  given,  we  need  the  appropriate  elastic  Green’s  function, 
known  for  only  a few  simple  shapes.  Since  we  know  the  image  tract ion 
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Pljn j but  not  \x^  on  S we  cannot  in  general  calculate  the  image  field 
in  the  body* 

The  change  in  volume  of  the  aolid  can  be  divided  Into  parts 
AV°°  and  AV*  arising  from  the  two  terms  in  (5).  tie  have  at  once 


AV 


oo 


f °° 

• I u . n 


dS  - 4*c  , 


(7) 


the  integral  being,  according  to  (1),  c times  the  solid  angle  sub- 
tended by  the  surface  at  the  singularity*  This  result  can  be  seen 

at  once  for  the  sphere* in*hole  model-  When  the  sphere  is  inserted, 

«© 

the  volume  of  the  hole  increases  by  4*c,  and  since  div^  1111  0 out* 
side  the  hole,  this  increase  is  transmitted  unchanged  to  S.  Again 
for  the  ideal  mathematical  singularity  for  which  (1)  holds  for  all 
£ we  have  strictly  not  div  U*  * 0 but  rather 

div  u°"  - • c V*  (l/'f  ) * )»  ($) 


with  a delta*function  of  expansion  at  ■ 0.  A formal  volume  in- 
tegration gives  (7)* 

Although  we  cannot  find  the  image  deformation  in  detail, 

we  can  find  AV*  with  the  help  of  the  rule  that  the  volume  change  of 
a body  subjected  to  a distribution  of  surface  traction  T_  per  unit 

area  is 


where  p^ 
For 


is  the  stress  produced  by  T and  K is  the  bulk  modulus* 


f (9,) 
j4.lV  u 
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since  :r  and  1x1  ela3tiic  equilibrium 

To  find  AV*  we  must  put  71  * ft.  /y\  . r - j>  . 71.* 

, oo  4 3*j3_’Jw 

with  pi . ' from  (4)„  But  the  operator  x^(  Q f o /L±  ) 

applied  to  H4.fe  * /^f  merely  multiplies  it  by  -2.  Hence 

the  integral  (9)  is  equal  to  the  integral  (7)  times  I#l/ 3K  and  eo 


.~CL 


4*  c 


1 ♦ C 


(10) 


Thus  the  total  volume  change  on  introducing  a center  of  dilatation 
of  strength  e into  any  homogeneous  isotropic  body  with  a stress-free 
surface  is 

AV  - AV~  ♦ tV1  - ^AV***  - 4sc  (11) 

where  the  constant 

V-3  J Ui2L  (12) 

1 ♦ <T 

is  about  1.5  for  metals  ( d~rw/l/3)  and  1.8  for  alkali  halides 
( <5“^  1/4) • This  result  has  already  been  given  by  Seitz2a. 


Frederick  Seitz,  Rev.  Hod.  Phys.  18,  384  (1946) 

It  is  also  implicit  in  the  known  result  for  the  interaction  energy 
of  a center  of  dilatation  and  an  external  hydrostatic  pressure2**: 
This  energy  is  the  product  of  AV  and  the  hydrostatic  pressure. 


2b 


See,  for  example,  J.  D.  Eshelby,  Phil.  Trans.  Roy.  Soc.  A 244.  87 
(1951) 
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If  In  (9*)  we  take  p^  to  be  any  state  of  purely  internal 
stress,  so  that  the  surface  integral  is  aero  we  arrive  at  the  known 
result  that  the  volume  of  a body  is  unaffected  by  internal  stresses 
provided  Hooke's  law  is  valid*  The  volume  expansion  we  have  calcu- 
lated is  not  in  disagreement  with  this*  Consider  the  hole- in-sphere 
model  of  section  2,  and  suppose  for  clearness  that  £>  is  negative* 
Then  the  unstrained  state  is  a body  with  a spherical  hole  of  radius 
rq  containing  a sphere  of  radius  r0(l  - (£,}  ) and  an  unoccupied 
volume  4* |£)  • 4*  jcj  ~j  - In  the  strained  state  the  surfaces 
of  the  hole  and  sphere  have  been  drawn  together  and  welded.  Slnoe 
the  total  volume  of  material  is  unchanged,  the  volume  enclosed  by 
the  outer  surface  muat  have  decreased  by  an  amount  equal  to  the 
volume  eliminated  between  sphere  and  hole.  This  simple  derivation 
does  not  supersede  our  more  elaborate  treatment  since  to  relate  the 
expansion  to  stress  effects  arising  from  the  imperfections  we  need 

oO  T 

to  know  something  about  and  separately*  If  there  are  depar- 
tures from  Hooke's  law,  the  volume  change  of  the  material  is  not 
quite  zero:  this  is  discussed  in  section  6. 


4*  The  Deformation  produced  by  a Large  Number  of  Singular  it  lea 

If  there  are  a number  of  singularities  in  the  body  we  have 


«» 


u 


(r)  *^cn 


X ~J& 

Z--**  1 


■3 


(13) 


with  summation  over  all  imperfections.  The  image  stress  p]j  is  that 
stress  which  has  no  singularities  within  3 and  on  S annuls  the  sur- 
face traction  calculated  from  (13).  u1  is  the  associated  displace- 


ment* 


If  the  defects  are  all  alike  we  have 


AVA  cf  /XX  (12* ) 

where  f is  the  atomic  fraction  of  defects  and  XX  is  the  volume 
per  atom  of  the  matrix. 

Of  special  interest  is  the  case  where  there  are  a large 
number  of  identical  centers  spread  through  the  volume , with  an 
approximately  uniform  density  of  n centers  per  unit  volume.  Here 
vre  can  also  say  something  about  the  change  of  shape  if  we  allow 
ourselves  to  replace  the  actual  distribution  by  a continuous  dis- 
tribution of  infinitesimal  centers  with  the  same  total  strength  per 
unit  volume.  (We  shall  try  to  justify  this  below}.  If  as  before 
the  body  is  embedded  in  an  Infinite  matrix  this  averaged  displace- 
ment is  clearly 

m 00 

u ■ cn 

at  any  rate  outside  S,  (It  is  not  at  once  clear  that  we  can  give 
any  meaning  to  the  average  displacement  in  a medium  riddled  with 
singularities.} 

The  comparison  already  made  with  thermal  expansion  sug- 
gests that  a body  uniformly  filled  with  infinitesimal  centers  of 
dilatation  will,  like  a uniformly  heated  solid,  expand  uniformly. 
The  following  is  a formal  proof.  We  first  show  that  the  displace- 
ment given  by  the  volume  integral  (14)  can  be  duplicated  outside  S 
by  a distribution  of  body  force  over  S of  uniform  magnitude 

p • K A 7/7  • 4#cnyK 

per  unit  area  and  directed  along  the  outward  normal  at  each  point. 

A point  force  F in  an  infinite  medium  produces  a displaces  t 
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a 

*w\ 


a L ♦ B £ (£>E> 
r 


b - i/i6^]i  (i  - <r  ) , 


A * (3  - 4<r  )B 


at  a vector  distance  r from  it.  Each  element  of  the  surfar.  * contri- 
butes  a force  dF  * p n d3  and  the  net  displacement  is 


Suppose  now  that  3 is  marked  out  in  an  infinite  medium, 
that  the  uniform  distribution  of  infinitesimal  centers  of  dilata- 
tion is  introduced  and  that  a layer  of  body  force  equal  and  opposite 
to  that  just  described  is  distributed!  over  3*  The  combined  effect 
of  these  last,  two  steps  gives  zero  displacement  outside  3 and  leaves 
3 its  original  chape  and  size.  If  the  unstrained  matrix  is  now  cut 
away,  the  elastic  state  of  the  body  bounded  by  3 is  unchanged,  but 
the  layer  of  body  force  becomes  a hydrostatic  pressure  p acting  on 
the  surface.  When  this  is  removed  S undergoes  the  uniform  expansion 

“i  ‘ -J*  n Xl  (15) 

o,-  Ci>  *r  d>*  t 

Though  u + u1  represents  a uniform  expansion,  u and  uA  separately 
do  not,  except  for  a spherical  body.  Thus  omission  of  the  image 
terms  in  addition  to  giving  a volume  change  lacking  the  factor  ^ 
{which  could  be  absorbed  in  the  usually  unknown  constant  c ) would 
also  give  0 non  •uniform  expansion. 


Y 
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Suppose,  for  example,  that  S is  an  ellipsoid  with  axes 
Then  u00  « - cn  grad  C j where 

* l/2  (D  - Ax^  - Bx^  - Cx?  ) (16) 

with  the  usual  notation  for  the  potential  of  an  ellipsoid.  Thus  if 
image  terms  are  neglected,  the  ellipsoid  becomes  another  ellipsoid 
with  semi-axes  * A s etc.,  where 


X 


A 

...  « cnA,  A & - cnB,  ~ — * cnC  . (17) 

tJl  J,  t 

Since  A ♦ B + C « 4*  the  volume  change  agrees  with  (7).  As  an  ex- 
ample, if  “ 2:1:1  then  A:B:C  « 2:1:1  very  nearly.  This  de- 

parture from  uniform  expansion,  if  it  existed,  would  be  hard  to  de- 
tect by  macroscopic  measurements,  but  as  we  shall  see  below,  there 
is  an  exactly  similar  effect  on  the  lattice  constants  determined  by 
x-ray  diffraction,  and  this  would  be  easily  observable. 

We  must  now  try  to  justify  replacing  the  sum  (13 ) by  the 
integral  (14).  Just  this  question  arises  when  we  replace  the  field 
of  a set  of  gravitating  or  electrified  particles  by  the  field  of  a 
continuous  body.  This  kind  of  problem  can  be  attacked  in  various 
ways,  none  completely  satisfying;  the  following  seems  a plausible 
line  for  the  elastic  case. 

The  observable  macroscopic  displacement  at  a point  of  a 
body  in  which  the  displacement  fluctuates  on  a microscopic  scale 
is  essentially  the  displacement  of  the  centroid  of  a region  large 
compared  with  the  scale  of  the  fluctuations.  It  is  therefore 
reasonable  to  define  the  gross  (or  macroscopic)  displacement  at  a 


12 


point  as  the  average  of  the  actual  (microscopic)  displacement  taken 
over  the  volume  of  a sphere  whose  radius  H is  large  enough  to  con- 
tain many  fluctuations , but  is  small  compared  with  the  dimensions 


of  the  body; 

Let  us  find  the  gross  displacemt  corresponding  to 
{14)  for  a point  inside  the  body.,  The  volume  integral  over  the 
sphere  E of  the  displacement  due  to  a singularity  at  r relative 
to  its  center  is  equal  to  the  attraction  on  a charge  c at  f' 

* s+~ 

produced  by  a distribution  of  charge  of  unit  density  filling  the 
sphere , i6e.  4/3  «cr  if  r 4 R and  4/3  *<3&r/r$  if  r > R*  Thus 


— f.Q 
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4.  y 

r3  *— 

"t<R 


£n 
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'OS. 


Jn 
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As  R increases  and  becomes  reasonably  large  compared  with  T] 
the  first  term  should  be  small,  being  proportional  to  tht  position 
vector  of  a large  number  of  points  taken  at  random  in  a sphere* 

At  the  same  time,  as  R increases  it  becomes  more  and  more  reason- 
able to  replace  the  second  term  by  an  integral  since  the  distance 
of  each  singularity  from  the  point  at  which  their  effects  are  summed 
is  large  compared  with  the  mean  distance  between  them.  V*e  thus 
reach  (14)  but  with  a sphere  of  radius  R about  the  point  of  observa- 
tion excluded  from  the  integration.  However,  by  a known  result  of 
potential  theory  (14 ) is  unaltered  by  this  omission  provided,  the 
sphere  lies  wholly  within  S„ 

Having  defined  the  gross  displacement  we  can  calculate  a 
gross  stress  from  it  with  the  aid  of  (3),  This  should  be  u.a  true 


macroscopic  stress,  i.e.  it  should  give  correctly  the  surface  trac- 
tion required  to  prevent  relative  motion  of  the  faces  of  a macro- 
scopic cut  made  in  the  material. 

We  can  obviously  generalize  these  results  to  a distribu- 
tion of  imperfections  in  which  is  a function  of  position.  The 
microscopic  u is  still  given  by  but  the  gross  displacement  is 


u 


C O 


n(r)  r 


dv 


(IB) 


and  in  the  averaging  process  we  must  have  a-l/3  4C  a « n-1/3 


|grad  (n_1/3 ) j . 

x 

The  gross  or  microscopic  Pjj  is  the  stress  produced  by 

OO 

surface  tractions  -PijUj  calculated  from  (13)  or  (1 4,  18),  These 
two  surface  tractions  will  differ  by  rapidly  fluctuating  quantities 
with  a ’wave-length 1 of  the  order  of  n-^/3  whose  effect  will  be  con- 
fined to  a surface  layer  of  the  same  order.  Thus  in  the  bulk  of  the 
material  the  image  quantities  are  efficiently  smoothed  by  1 f » 

Venant’s  principle  whether  we  smooth  the  infinity  quantities  or  not  } 
and  we  need  not  distinguish  between  their  gross  and  microscopic  values^ 


5 • The  Effect  of  Point  Singularities  on  the  X-ray  Diffraction  Pattern 
Miller  and  Bussell  originally  suggested  that  there  Is  a 


3 

P.  H„  Miller,  Jr„,  and  B.  R.  Russel,  J.  Appl.  Phys.  ,23,  II63  (1952), 


difference  between  the  geometrical  volume  change  of  a crystal  contain- 
ing imperfections  and  the  volume  change  deduced  from  the  change  of 


x-ray  lattice  constant.  They  calculate  that  for  a uniform  distri- 
bution of  point  imperfections  the  latter  should  be  about  tv/ice  the 
former, 4 Huang5  found  that  they  were  equal.  These  authors  omit 

^ They  have  since  withdrawn  this  factor,  P,  H,  Miller,  Jr«  and  B„ 

R«  Russel,  J*  Appl«  Phys. 

5 K.  Huang,  Proc„  Roy.  Soc.  A190,  10 2 (1947)o 

image  effects.  Miller  and  Russel  calculate  the  image  effect  for  a 
single  imperfection  at  the  center  of  a sphere.  They  take  » collec- 
tion of  such  spheres  with  the  interstices  filled  up  to  represent  a 
body  containing  many  imperfections  and  argue  that  the  image  effect 
is  negligible  because  the  surface  traction  across  these  spheres  is 
by  no  means  zero,  As  we  have  seen,  the  effect  of  the  image  term 
on  the  volume  change  can  be  taken  into  account  quite  generally, 
and  we  shall  fine  that  tlie  3ame  is  true  for  its  effect  on  the  x-ray 
lattice  constant* 

For  an  undisturbed  crystal  the  scattering  power  plotted 
in  reciprocal  space  will  be  a set  of  patches  with  their  maxima  at 

the  points  of  the  reciprocal  lattice*  If  the  crystal  is  distorted 
by  external  forces  or  internal  imperfections  the  patches  will  be 

;:  "'jac*Hi  u or::.*M  n'<-  y r .*fard  the  new  maxima  as  defining  a 

m - a ,;.o-  a.  tt  and  wo  car.  >.  ht:c  v-pe,;.  of  the  deformation  of 
the  re-,  ipr  ' ai  la . t ice  corresponding  to  a deformation  of  the  crystal, 
let  the  position  cf  a point  of  ? he  crystal  or  .reciprocal 
lattice  before  and  after  distortion  be 
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Miller  and  Russel  derive  the  following  relation  valid  within  a re- 
gion around  the  origin  of  reciprocal  space  which  decreases  as  the 
degree  of  crystal  distortion  increases: 

Ahi  Btj  (j  - 1,2,1)  (19) 

where 

Aij  “ Li  BiJ  " Lj» 


implying  summation  over  all  points  of  the  crystal  lattice*  We 
have  simplified  equations  (1,  2,  3)  in  reference  3,  appendix  1,  by 
taking  the  origin  at  the  lattice  point  nearest  the  center  of  gravity 
of  the  crystal:  then  the  sums  will  differ  from  zero  by  inte- 

gers small  compared  with  the  number  of  atoms  in  the  crystal  and  may 
be  neglected.  Equation  (19)  has  the  solution 


Ah, 


* hk  Cki» 


where  A . is  the  matrix  reciprocal  to  A 

*3  ij 


°ki  * Bkj  Aji  (20) 

C is  not  in  general  sym- 
metric even  if  B^j  is.  Thus  (20)  is  an  infinitesimal  affine  trans- 
formation in  reciprocal  space.  In  other  words  the  reciprocal  lattice 
points  arc  shifted  as  if  they  were  embedded  in  am  imaginary  elastic 
continuum  which  undergoes  a deformation  in  which  all  the  strain  and 
rotation  components 


•ij  * | ( 0 tol/B  hj  * 9a h}/9  hi)  - (Cij  ♦ Cji ) 

W - k ( 9 hj  - - -ifCj,  ..  G") 

2 2 
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are  in  general  different  from  zero  and  constant  from  small 
This  implies,  for  example,  chat  from  x-ray  measurements  of  the 
position  of  low-order  spots  alone,  a distorted  cubic  crystal  could 
not  be  distinguished  from  an  unstrained  monoclinic  crystal  slightly 
misoriented* 

Following  Huang  we  may  replace  the  displacement  expressed 
as  a sum  of  the  effects  of  the  separate  imperfections  by  a*,  integral 
if  the  distribution  is  statistically  uniform*  This  corresponds  to 
the  averaging  process  of  the  previous  section,  and  we  may  put  aaLi  • u^. 
The  summation  may  also  be  replaced  by  an  integral  and  we  have 

x±  x<  dv,  Bij  " ^5  j *i  Uj  dv  {21) 

But  we  know  that  when  the  image  terms  are  included  has  the  form 
(15),  so  that  A, ^ * Jtft  n c B - . and  (20)  has  the  solution 

Ahi  Aho  Ah-j  r 

— — “ m — »*■  * •"  **  *ncy 

hl  h2  h3  3 * 

(The  A,,  are  the  products  of  inertia  of  the  crystal  and  dat  (A. 4) 
ij 

could  only  vanish  if  the  moment  of  inertia  of  the  crystal  about  some 
axis  were  zero, ) Thus  the  reciprocal  lattice  contracts  uniformly 
in  the  way  we  should  expect  for  a crystal  which  had  suffered  the 
Tin  if  orm  expan  si  on  ( 15 ) - 

Thi3  simple  result  would  be  destroyed  if  the  image  terms 
were  omitted.  We  should  have 
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where 

f * X dy/r 

Is  the  potential  of  the  specimen  if  filled  with  charge  of  unit  density. 
For  an  ellipsoidal  crystal  CP  is  a quadratic  function  of  the  coordi- 

I ' * /n 

nates  (eqn„  {16}),  and  it  follows  that  -Ahj/h^,  *^3 A3  are 

equal  respectively  to  the  three  quantities  (17  )•  In  this  case  the 
(incorrect)  deformation  of  the  reciprocal  lattice  agrees  with  the 
(incorrect)  macroscopic  deformation  with  image  effects  omitted.  For 
a body  of  arbitrary  shape  this  is  no  longer  true,  since  there  is  no 
simple  relation  between  Ajj  and  B-jj.  X-ray  measurements  would  in- 
dicate that  the  cubic  unit  cell  had  become  monoclinic,  and  its  change 
of  volume  would  bear  no  simple  relation  to  the  change  in  volume  of 
the  crystals  calculated  with  or  without  image  effects.  These  com- 
plications are  implicit  in  Huang's  calculations.  Each  of  the  two 
terms  his  expression  (I)  is,  in  electrostatic  language,  the  sum  of 
the  potenti alt'  of  a set  of  parallel  dipoles.  Such  a sum  taken  over 
a large  finite  region  in  an  infinite  distribution  of  dipoles  is 
ultimately  independent  of  the  size  of  the  region  but  depends  on  its 
shape  and  the  relative  position  of  the  point  of  observation  within 
it  . Inclusion  of  the  image  terms  adds  a layer  of  surface  charge 
which  just  annuls  these  variations.  Thus  Huang's  (I)  does  not  re- 
rrpf.ent.  a uniform  expansion  > .-<*  t for  a sphere,  for  the  latter  the 
ex’S-r.rior:  locks  factor  ''  ' y r :in>r  c-  .,<  > t ; . ' a the 

- ■ ; ; ; .1  :,f  .<.u  w.,-  • • 

*.  ;*•»>•;  ’•<  . m iai  t-ice  : » * ,y  -e**  . - *•  ■ r.  • . • 

■ I ) a • III)  1.  !•-.  t hi-  > 1 ‘jy  c J , :i.‘  i--  r !s 
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eventually  eliminated  with  the  help  of  his  equation  (IS)  in  which 
also  c mu3t  be  replaced  by  the  final  result  is  unaffected. 

Miller  and  Russel  found  by  calculating  the  5. 
numerically  that  for  a single  singularity  at  the  center  of  a sphere 
the  expansion  deduced  from  x-rays  was  2*5  times  the  geometrical  ex- 
pansion* They  deduce  from  this  that  there  should  be  a similar 
difference  {with  the  factor  reduced  from  2*5  to  about  2)  when  there 
are  a large  number  of  imper feet  ions  distributed  at  random* 

The  previous  discussion  indicates  that  this  is  incor- 
rect, but  does  not  show  exactly  how  the  discrepancy  arises * This 
discrepancy  exists  whether  image  terms  are  introduced  or  not (though 
introducing  them  would  reduce  the  factor  from  2 to  (1  ♦ ^ )/  ^ ) 

and  for  simplicity  we  omit  them  in  the  following  discussion. 

Let  us  find  the  effect  of  a single  imperfection  at  a 
point  7*  relative  to  the  center  of  a spherical  crystal  of  radius 
R.  We  can  save  calculations  by  finding  the  dilatation  of  the  re- 
ciprocal lattice  and  not  Its  detailed  distortion.  Putting 
Ut  • c(x^  • £ ^)/  | x - | in  {21,'  we  find  after  some  man- 

ipulation 
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sphere  of  unit  density  and  has  the  value  « (}H 
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If  we  exclude  from  the  integration  a sphere  of  radius  rQ  about 
a term  »r0  /R  must  be  introduced  inside  the  bracket.  Since  r0 
can  presumably  be  taken  of  atomic  dimensions  we  shall  omit  it. 

The  reason  for  the  discrepancy  now  becomes  clear.  The 
expression  (22)  is  greater  or  less  than  the  geometrical  AY  ac- 
cording as  4*  is  greater  or  less  than  ^(3/5 )&  * Q.77R.  For 

* 0 it  is  2.5  times  the  geometrical  value,  in  agreement  with 
Miller  and  Russel.  If  there  are  a large  number  of  defects  scat- 
tered at  random  throughout  the  sphere,  the  average  contribution 

06  i i- 

to  AV  per  imperfection  is  given  by  replacing  l in  { 22 ) by  its 
mean  value  for  all  points  of  the  sphere.. 


4# 


d.yV* 
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r dr 


ihis  gives  a mean  x-ray  A'/  of  4*e,  in  agreement  with  the  geome- 
trical value  {7}«  There  seems  to  be  nothing  objectionable  it:  the 

replacement  of  sums  by  integrals,  and  presumably  Miller  and  Russel fs 
summation  procedure  would  give  the  geometrical  value  if  carried  out 
for  a reasonable  number  of  imperfections  scattered  throughout  a 
sphere. 


As  a corollary  we  see  that  with  a given  number  of  defects 
the  x-ray  expansion  depends  strongly  on  their  (non-uniform)  distri- 
bution, while  the  geometrical  expansion  is,  of  course,  independent 
of  the  distribution.  Thus,  for  a sphere  the  ratio  of  the  x-ray 

oO  aft 

AY  to  the  geometrical  AY  can  be  given  any  value  between  2.5  and 
0 by  sweeping  the  defects  towards  the  center  or  surface.  When 
the  two  expansions  are  corrected  for  the  image  terms  the  ratio  lies 

i 8 


between  (1.5  ♦ 
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and  ( 
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6«  Discussion 

The  equality  of  the  x-ray  and  geometrical  expansions  of  a 
crystal  uniformly  filled  with  point  imperfections  has  a bearing  on 
the  question  whether  they  are  of  Schottky  or  Frenkel  typc3f  Revised 
calculations  by  Miller  and  Ruseel^  gave  a difference  of  nci  more 
than  10  to  20  per  cent  difference  between  the  two  expauaicrs  in 
place  of  their  original  factor  of  two,.  The  residual  difference 
raay  be  partly  due  to  the  neglect  of  image  terms. 

Turning  to  the  factor  ^ introduced  into  the  expansion 
by  the  image  terms  it  looks  at  first  sight  as  if  one  could  rede- 
fine the  strength  of  an  imperfection  as  c and  then  forget  about 
the  image  effects.  This  is  not  so.  The  contributions  to  the  ex- 
pan  si  on  from  u and  vlx  are  in  the  ratio  1:  [y  - 1).  But  -;he 

*kv-  -w-  “ 

ratio  of  their  contributions  to  any  other  phenomenon  will  in  gener- 
a!  be  quite  different,  since  u is  a rapidly  fluctuating  function 
of  position  associated  with  a pure  shear  whereas  is  smooth  ex- 
cept near  the  surface  with  an  associated  stress  which  is  chiefly  a 
hydrostatic  pressure  in  a body  of  reasonable  shape,  c is  thu3  the 
basic  constant  which  we  might  hope  to  determine  from  measurements 
of  one  effect,  and  apply  in  calculating  the  value  of  another.  The 
following  examples  illustrate  this. 

Rabarro^  has  discussed  the  hardness  of  dilute  alloys. 
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F.  a.  N.  Nabarro.  hroc.  fhys.  Soc.  (London)  £6  , 669  (1946), 


Oo 


The  hardening  is  la.'  entire] y to  p^  , the  slowly  varying 

pu  playing  scarcely  any  part.  He  eliminates  c with  the  h-.p  of 
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the  observed  change  of  lattice  constant  with  composition  but  neglects 
the  image  effect.  This  can  be  put  right  by  dividing  his  theoretical 
constant  oo*  0,05  by  ^ before  comparing  it  with  the  values  deduced 
from  experiment.  The  spread  of  the  latter  is  too  great  to  say  if 
this  gives  any  improvement, 

Dexter?  has  calculated  the  effects  of  u and  u*on  the 

■ ■ I»iwi’«'iwii  ■!  I I I *«■»■  *■■■■»  ■■>■  ■— t — ■ lA'nMwKxnna'  imuwinwi'i  I I — mwi  m wi  ■winn^.,i] 

7 D.  L.  Dexter,  Phys.  Rev.  8 £,  768,  1952, 

electrical  conductivity  of  metals.  The  former  gives  a tecaperature 
independent  and  the  latter  a-  temperature  dependent  contribution. 

They  are  of  the  same  order  at  room-temperature, 

Overhaus er^  has  considered  the  modification  of  ratio  of 


g 

A,  W.  Overhauser,  Phys.  Rev.  Q0 , 3C)3  ('993s) 

diffusion  by  the  strain  due  to  point  imperfections.  The  effects 
arising  from  u and  u!  depend  differently  on  temperature,  and  the 
image  term  may  become  important  at  high  temperatures* 

Zen er^  has  calculated  the  change  of  shear  modulus  caused 

^ C.  Zener,  Acta  Cryst.  2,  163  {1949). 

by  the  u field  of  point  imperfections.  Evidently  tha  dilatation 
due  to  the  image  field  should  give  rise  to  an  analogous  change  of 
bulk  modulus.  (As  i/e  have  seen  u*°  produces  no  dilatation  of  the 

*+*S.*M  *■ 

material  between  the  imperfections, ) In  fact  we  should  have  h change 

AK  _ d Iof.  K 6V1  _ d log  K / - 1 du 

K d log  V V d log  V * y 7 <2^ 


proportional  to  the  relative  change  d {^/\  of  lattice  constant  but 
of  opposite  sign,  since  d log  K'/ k log  V is  negative.,  According  to 
a simple  theory^-®  verified  quite  well  by  experiment  ~d  log  K/d  log  V 


J*  C.  Slater,  Phys*  ilev.  J&;1  744  (1940). 

is  equal  to  2 P + 4/3,  where  "P  is  Gruneisen’s  constant,,  Zener’s 
on  the  other  hand  depends  an  the  shear  strain  energy  associated 
with  an  imperfection  and  so  is  proportional  to 

By  an  atomic  calculation  Dienesll  has  found  that  one  per 
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G.  J*  Dienes,  Phys,  Rev.  66,  225  (1952)* 


cent  of  interstitials  increase  the  bulk  modulus  of  copper  by  6*6 
per  cent  and  of  sodium  by  1„9  per  cent,  His  calculations  deal  in 
principle  with  vdiat  happens  in  the  immediate  neighborhood  of  the 
imperfection,,  so  it  is  reasonable  to  suppose  that  they  take  account 

of  u and  that  the  effects  of  ux  should  be  added*  The  results 

- 

vcnA 

of  Dienes  and  Huntington  and  Seltz^*2  suggest  that  we  may  take 
12 

"H.  B,  Huntington  and  F,  Seitz,  Phys.  Rev,  <j>l,  315  (1942 }„ 

Ws/JTjL  t0  be  about  0.6  and  0.9  respectively  for  an  interstitial 
in  copper  or  sodium.  By  (12*)  and  (23)  the  corresponding  percen- 
tage changes  for  one  per  cent  interstitials  are;  lattice  constant 
0,9  (Cu),  1.4  (Na),  bulk  modulus  - 4.6  (Cu),  - 6,9  (Na).  Our 
calculation  cannot  claim  any  accuracy,  but  it  suggests  that  the 
image  term  is  important,  perhaps  even  large  enough  to  change  the 
3igp  oi*  '.••K/K, 


V/ 


The  change  of  shear  modulus  fvund  by  Zener’s  argument 
depends  on  non-linear  behavior  in  the  is  gion  close  to  the  .Imperfec- 
tion (see  belovi)  and  it  is  therefore  in  principle  already  included 
in  a calculation  such  as  Dienes*  which  co>' aiders  the  balance,  of 
interatomic  forces  near  the  imperfection* 

It  was  pointed  out  at  the  end  of  « action  3 that  for  a 
purely  elastic  model  the  actual  voliune  of  k » terial  was  unaltered* 
Zener 13  has  shown  how  in  such  a case  the  departure  from  Ik  oka's 


C*  Zener,  Trans-  Am-  Inst-  Mining  Met,,  Eng.  1^7,  361  (1942 h 


law  gives  a volume  change,  and  we  ought  to  consider  what  effect 
this  will  have  on  our  calculations*  This  change  has  the  value 
4V  - ~(1  + d log  K/d  log  V)  rfd  A 

~(1  ♦ d log  /d  log  V)  kjA 


where  and  are  the  total  dilatations i and  shear  strain  ener- 
gies associated  with  the  internal  stress*  The  first  term  can  be 

T 

considered  as  a small  correction  to  corresponding  to  i sing  in 

(9)  the  value  of  K appropriate  to  the  lattice  as  expanded..  The 

shear  strain  associated  with  the  sphere- in -hole  model  of  section 
x 

2 is  3*  m-  c /r0--  (of  ref*  3).  If  we  take  the  case  of  interstitials 
In  a face-centered  cubic  material  and  put  r0  equal  to  the  distance 
from  the  interstitial  to  the  nearest  face-centered  atom  m have 


?„ 


Th  is 
tv  for 


AV'A  * - (1  * d log  / d log  V)  ( j^/Y.)  (c/J  .1  > 

i a atxujfc  one  tenth  of  (12* ) for  int  erstitials  in  copj  er  if 
■e  v?  rake  c./_0-  -05*  However,  even  i.f  the  r-at  rk 
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cuito  l-irgfe  our  previous  calculations  would  l-c-t  be  up -set . 
shear  strain  energy  jwc  /r0^  associated  with  an  imperf:-.  .bio» 
nearly  ninety  percent  is  contained  within  a sphere  of  relive 
Thus  the  correction  coin:%s  entirely  from  an  additional  oxpaagicn  con* 

f ; 4 . ■ ' ) 

fined  to  the  immediate  neighborhood  of  the  imperfection  . 1 -‘‘his’  ex- 
tra volume  change  is  transmitted  to  the  surface  by  the  displacement 
which  must  have  the  form  (1)  as  soon  as  the  elastic  region  is 

tw- 

reached-.  Thus  the  effect  we  are  considering  merely  *renorsiaijses 1 
the  constant  c in  (1)  and  it  is  correct  to  use  the  modified  value 
not  only  in  (12*)  but  also  in  discussing  any  other  phenomenon  in 
which  only  strains  in  the  elastic  region* 

Point  imperfections  may  expand  the  lattice  by  mechanisms 
other  than  misfit,,  Take  fox’  example  an  alkali  halide  crystal  with 
doubly  charged  foreign  cations  in  some  positions  together  with  an 
equal  number  of  vacant  anion  sites*  The  continuum  model  is  a di- 
electric with  positive  and  negative  charges  embedded  in  it,.  Round 
each  there  is  an  elastic  stress  field  due  to  eiectroatr  lotion* 

The  dilatation  is  large  near  a change  and  falls  off  rapidly  with 
distance,,  The  stress  field  is  essentially  that  of  a center  of 
dilatations  A rough  calculation  suggests  that  its  strength  is  of 
the  same  order  as  for  a reasonable  degree  of  misfit..  It  is  of  the 
same  sign  for  both  positive  and  negative  charges*  Again,  if  atoms 
of  different  valency  are  introduced  aubstitutionally  into  a metal 
the  resulting  change  in  the  number  of  electrons  in  the  conduction 
band  will  produce  a change  of  lattice  constant  quite  apart  from  any 
misfit  effects* 


It  is  not  the  aim  of  this  paper  to  advocate  the  elastic 
approach,  but  to  shotv  that  when  it  is  used,  neglect  of  the  pre~ 
sence  of  the  free  surface  may  lead  to  qualitative  untidinesses 
(such  as  non-uniform  deformation  by  a uniform  distribution  of  de- 
fects) and  quantitative  errors  which  may  be  important Q 


